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ABSTRACT 

We present the discovery of 17 low mass white dwarfs (WDs) in short-period P < 1 day binaries. 
Our sample includes four objects with remarkable \ogg — 5 surface gravities and orbital solutions 
that require them to be double degenerate binaries. All of the lowest surface gravity WDs have metal 
lines in their spectra implying long gravitational settling times or on-going accretion. Notably, six 
of the WDs in our sample have binary merger times <10 Gyr. Four have >0.9 Mq companions. If 
the companions are massive WDs, these four binaries will evolve into stable mass transfer AM CVn 
systems and possibly explode as underluminous supernovae. If the companions are neutron stars, then 
these may be milli-second pulsar binaries. These discoveries increase the number of detached, double 
degenerate binaries in the ELM Survey to 54; 31 of these binaries will merge within a Hubble time. 
Subject headings: binaries: close — Galaxy: stellar content — Stars: individual: SDSS J0751-0141, 
SDSS J0811-f0225 — Stars: neutron — white dwarfs 



1. INTRODUCTION 

Extremely low mass (ELM) WDs, degenerate objects 
with \ogg < 7 (cm s~^) surface gravity or <0.3 Mq mass, 
are t he product of com mon envelope binary evolution 
(e.g. IMarshet aLlll995[ ). ELM WDs are thus the sign- 
posts of the type of binaries that are strong gravitational 
wave sources and possible supernovae progenitors. The 
goal of the ELM Survey is to discover and characterize 
the population of ELM WDs in the Milky Way. 

Previous ELM Survey papers have reported the dis- 
covery of 40 WDs spanning 0.16 Mq to 0.49 Mq found 
in the Hypervelocity Star (HVS) Survey, the Sloan Dig- 
ital Sky Survey (SDSS), and in our own targeted sur- 
vev (IBrown et al.i[2010ll20TIbll2012bl:IKilic et al.l[2010R 
I2011al l2012(r We refer to this full sample of WDs as 
the ELM Survey sample, but reserve the term "ELM 
WD" for those objects with log 5 < 7. All of our 
WDs are found in short-period, detached binaries, 60% 
of which have merger times <10 Gyr. Three notable 
systems are detached bina ries with <40 min orbital 
periods (|Kilic et al.l I2011bl l3). The eclipsing system 
J065 1+2844 is the second-strongest gravitational wave 
source in the mHz range (Brown et al. 201 la) . We mea- 
sured it s period change in on e year with optical eclipse 
timing (jHermes et al.l [2012bD . Other results from the 
ELM Survey include the first tidally distorted WDs 
(jKilic et al.ll2011cH Hermes et al."2012a') and the first pul- 
sating helium-core WDs (Hermes et al. 2012c, 2013). 

Here we present the discovery of 17 new WD bina- 
ries identified fr om spectra previously obtaine d for t he 
HVS Survey of IBrown et a l. (2005, ' 2006a, b], |2007a||3, 
I2009ll2012al) . iKilic et al.l " (J2b07a) analyzed the visually- 
identified WDs in the original dataset and discovered one 
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Based on observations obtained at the MMT Observatory, a 
joint facility of the Smithsonian Institution and the University 
of Arizona. 



ELM WD binary (|Kilic et al.ll2007bD . This approach 
failed to identify the lowest surface gravity WDs. We 
now fit stellar atmosphere models to the entire collection 
of spectra not previously identified as WDs, and acquire 
follow-up spectroscopy of new ELM WD candidates. The 
result of this effort is that we find low surface gravity ob- 
jects that might not be considered WDs if not for their 
observed orbital motion. 

We chose to call objects with 5 < logg < 7 "ELM 
WDs" because these objects occupy a unique region of 
surface gravity/effective temperature space that overlaps 
the terminal WD cooling branch for ELM WDs (see Fig- 
ure[T]). This region is well separated from both hydrogen- 
burning main sequence tracks and helium-burning hor- 
izontal branch tracks. Our objects are systematically 
^^10,000 K too cool, give n their surface gravi ties, to be 
helium burning sdB stars. iKaplan et al.l (|2013l ) refer to a 
similar type of object as a "proto-WD." The lowest grav- 
ity objects in our sample may indeed have hydrogen shell 
burning and thus are, properly speaking, proto-WDs, 
but here we address our sample of low gravity objects 
as ELM WDs. A variety of observations demonstrate 
that degenerate helium-core WDs exist at our observed 
temperatures and surface gravities. ELM WD compan- 
ions to milli-second pulsars are direc tly observed (e.g. 
iBassa et al.l 120061 : iCocozza et al.l 12006 ) at the tempera- 
tures and gravities targeted by the ELM Survey. The 
measured radius of the log 5 = 6.67 ± 0.04 object in 
J0651-F2844, a 0.0371 ± 0.0012 i^p, star, demonstrates it 
is a de g enerate WD (IBrown et al.l I2011cl : IHermes et al.l 
I2012b|) . iVan Grootel et al.l (|20T3f ) account for the un- 
usually long pulsation periods of the ELM WD pulsators 
with low mass WD models. Hence, calling our low sur- 
face gravity objects ELM WDs is appropriate. 

Interestingly, we find ELM WDs in binaries with >0.9 
Mq companions and rapid merger times. When these de- 
tached binaries begin mass transfer. Marsh et al. (200^ 
show that the extreme mass ratios will lead to sta- 
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ble mass transfer. For the case of massive WD accre- 
tors, theorists predict large hehum flashes that may ig- 
nite thermonuclear transients du bbed ".la" supernovae 
(jBildsten et al.ll2007l : IShen fc Bil dstcn 2009) or may det- 
onate th e surface helium-layer and t he massive WD 
(INomotd[T982: Wooslev fc Weaver 1991 ISim et al.l[20Tl 
and produce an underluminous supernova. However, the 
final outcome of such mergers is un certain and they may 
not trigger supernovae explosions ()Dan et al.ll2012f ). If 
the companions are instead neutron stars, this would be 
the first time a milli-second pulsar is identified through 
its low-mass WD companion. Such systems allow mea- 
surement of the binary mass ratio and the neutron star 
mass through a combination of the pulsar orbit obtained 
from radio timing and the WD orbit obtained from op- 
tical radial velocity observations. 

The importance of identifying the ELM WDs in the 
HVS Survey is that the HVS Survey is a well-defined 
and a nearly 100% complete spectroscopic survey. With 
a complete sample of ELM WDs we can measure the 
space density, period distribution, and merger rate of 
ELM WDs, and link ELM WD merger products to pop- 
ulations of AM CVn stars, R CrB stars, and possibly un- 
derluminous supernovae. In a stellar evolution context, 
our ELM WD survey complements studies of WD bina- 
ries with main sequence companions (e.g. iPvrzas et al.l 
l2012t iRebassa-Man sergas et al."2012i) and with sdB sta r 
companions (e.g. [Gcicr ct al. 2012; S ilvotti et al.l 120121 ). 
sdB stars are TpM > 25, 000 K helium-burning precursors 
to WDs (jHeberi 12009 ). Our ELM WDs, on the other 
hand, have T^eK 20,000 K. 

We organize this paper as follows. In Section 2 we 
discuss our observations and data analysis. In Section 
3 we present the orbital solutions for 17 new ELM WD 
binaries. In Section 4 we discuss the properties of the 
ELM WD sample. We conclude in Section 5. 

2. DATA AND ANALYSIS 
2.1. Target Selection 

The HVS Survey is a targeted spectroscopic survey of 
15 < go < 20 stars with the colors of ~3 Mq main se- 
quence stars, stars which should not exist at faint magni- 
tudes in the halo unless t hey were ejected ther e. The tar- 
get selection is detailed in lBrown et al.l ()2012aD and spans 
-0.4 < [g- r)o < -0.25, 0.4 < [u - 5)0 < 1.07. This 
color selection fortuitously targets WDs in the approxi- 
mate range 10, 000 < T^s < 20, 000 K and logg < 7.5. 

We separate ELM WDs from the other stars in the 
HVS Survey using stellar atmosphere model fits to 
the single-epoch spectra. Our initial set of fits uses 
an upgraded version of th e ferre code described by 
lAllende Prieto et al.l (|2006( ) and synthetic DA WD pure 
hydrogen spectra kindly provided by D. Koester. The 
grid of WD model atmospheres covers effective tempera- 
tures from 6000 K to 30,000 K in steps of 500 K to 2000 
K, and surface gravities from log 5 =5.0 to 9.0 in steps of 
0.25 dex. The model atmospheres are calculated assum- 
ing local thermodynamic equilibri um an d include both 
convective and radiative transport (jKoest cr 2008). 

W e fit 2,000 spectra, niostly f ro m the o ri ginal HVS Sur- 
vey (IBrown et al.ll2005l l2006allbl. I2007allb| . 120091) but also 
some newly acquired data (jBrown et al1l2012a[) . These 
2,000 spectra were not previously analyzed for the ELM 
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Fig. 1. — Surface gravity vs. effective temperature of the 57 WD 
candidates (blue diamonds) and the previously published ELM 
Survey WDs (red circles). The 17 newly discovered binaries are 
plotted with errorbars and parameters obtained from Gianninas 
stellar atmosphere models; the parameters of the other WDs were 
obtained from Koester models. For reference, we plot theoretical 
tracks for .16-0.45 A/0 hydrogen atmosphere WDs (green lines, 
IPanei et al.l 2007), main sequence tracks for 2, 3 and 6 Mq ' 



(cyan lines. [Girardi et aEII200J ). and horizontal br anch tracks for 
0.488 , 0.495, and 0.500 Mq stars (magenta lines, IDorman et al.l 
119931) . Zero-age main sequence (ZAMS) and zero-age horizontal 
branch (ZAHB) isochrones are drawn with thick lines, a s is the ho- 
mogen ous helium-burning main sequence (He-MS) from lPaczvhskil 

mm . 

Survey because the spectra were not previously identi- 
fied as WDs. We fit both flux-calibrated spectra (for im- 
proved Toff constraints) as well as continuum-corrected 
B aimer line profiles (insensitive to reddening and flux 
calibration errors). For the continuum-corrected Balmer 
line profiles, we normalize the spectra by fitting a low- 
order polynomial to the regions between the Balmer 
lines. We adopt the parameters from the fiux-calibrated 
spectra, except in cases where the spectra were obtained 
in non-photometric conditions. The uncertainties in our 
single-epoch measurements are typically ±500 K in Tcb 
and ±0.1 dcx in logy. From these fits we identify 57 low 
mass WD candidates. 

2.2. New Spectroscopic Observations 

We obtain follow-up spectra for each of the candidate 
low mass WDs to improve stellar atmosphere parameters 
and to search for velocity variability. 

Observations were obtained over the course of seven 
observing runs at the 6.5m MMT telescope between 
March 2011 and Febru ary 2013. W e used the Blue Chan- 
nel spectrograph (.Schmidt et al.| [l989') with the 832 line 
mm~^ grating, which provides a wavelength coverage 
3650 A to 4500 A and a spectral resolution of 1.0 A. 
All observations were paired with a comparison lamp ex- 
posure, and were fiux-calibrated using blue spectropho- 
tometric standards (Masscv ct al. 198^. The extracted 
spectra typically have a signal-to-noise (S/N) of 7 per 
pixel in the continuum and a 14 km s~^ radial velocity 
error. 

We obtained additional spectroscopy for g < 17 mag 
ELM WD candidates in queue scheduled time at the 
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Fig. 2. — Gianninas model fits (smooth red fines) overplotted 
on the composite observed spectra (black lines) for the 17 WD 
binaries. The Balmer lines are arranged from H12 (top) to H7 
(bottom); left panel plots the K < 200 km s~^ binaries, right 
panel plots the K > 200 km s^^ binaries. 



1.5m FLWO telescope . We used the FAST spectrograph 
(jFabricant et al.lll998D with the 600 line mm"^ grating 
and a 1.5" slit, providing wavelength coverage 3500 A 
to 5500 A and a spectral resolution of 1.7 A. All obser- 
vations were paired with a comparison lamp exposure, 
and were flux-calibrated using blue spectrophotometric 
standards. The extracted spectra typically have a S/N 
of 10 per pixel in the continuum and a 18 km s^^ radial 
velocity error. 

2.3. ELM WD Identifications 

Our follow-up observations provide improved stellar at- 
mosphere constraints, from which we determine that 24 
(42%) of the candidates are probable ELM WDs with 
5 < \ogg < 7 (see Figure [1]). The remaining candi- 
dates have either log 5 > 7 and are normal DA WDs, or 
log 5 < 5 and are presumably halo blue horizontal branch 
or blue straggler stars. 

Figure [1] plots the distribution of T^s vs. log 5 for all 
of the WDs with logg > 5. Previously published ELM 
Survey stars are marked with red circles. The observed 
WDs overlap tracks based on Pan el et al.l (pOOT) mod- 
els for He-core WDs with hydrogen shell burning ((green 
lines) iKilic et aLll2010lJ ). but do not overlap main se- 
quence nor horizontal branch evolut i onary tracks. For 
reference, we also plot iGirardi et al.l (|2002L 120041 ) solar 
metallicity main sequence tracks for 2, 3, and 6 Mq 
stars (cyan lines). A helium-burning horizontal branch 
star can h ave a surface grav ity similar to a degenerate He 
WD fe.g.. lHeber et al.|[2003f ). but only at a systematically 
higher effective temperature than that targeted by the 
ELM Survey. This is illustrated bv lDorman et al.l ()1993l ) 
[Fe/H]=-1.48 horizontal branch tracks for 0.488, 0.495, 
and 0.500 Mq stars (magenta lines), as well as the ho- 
mogen ous helium-burning main sequence from iPacz vhskil 
()1971|) . The zero-age main sequence and helium-burning 
horizontal branch isochrones (thick lines) mark the lim- 
its. 

Seventeen of the WDs show significant velocity vari- 
ability, including twelve of the newly identified ELM 
WDs. The other ELM WDs have insufficient coverage for 
detecting (or ruling out) velocity variability. We focus 
the remainder of this paper on the 17 well-constrained 
systems. 



2.4. Improved Atmosphere Parameters 

Our sample of 17 well-constrained systems contains 
four logg ~ 5 objects, objects at the edge of our stel- 
lar atmosphere grid that might not be considered de- 
generate WDs if not for their obs erved or b ital m otion. 
We also fit these stars to grids of iKuruc^ ()1993l ) mod- 
els spanning —4 < [Fe/H ] < 0.5, 1 < logg < 5 (see 
lAllende Prieto et"aLl 120081 ) but found very poor solu- 
tions. This motivated us to perform a new se t of stellar 
atmosphere fits using iGianninas et al.l ()2011[ ) hydrogen 
model atmospheres. The Giannin as models employ im- 
prove d Stark broadening profiles (jTremblav &: BergeronI 
I200 9D. and the ML2/a = 0.8 prescription of the mixing- 
length theory for models wher e convective energy trans- 
port is important t Tremblav et al.ll201Ql . We calculate 
atmosphere models for low surface gravities, and use a 
grid that covers Tes from 4000 K to 30,000 K in steps 
ranging from 250 to 5000 K and logg from 5.0 to 8.0 in 
steps of 0.25 dex. 

Our stellar atmosphere fits use t he so-called spectro- 
scopic technique described in iGianninas et al.l (|2011). 
One d ifference between our work and IGianninas et al.1 
()2011|) is that we fit higher-order Balmer lines, up to 
and including H12, observed in the low surface gravity 
ELM WDs. The higher-order Balmer lines are sensitive 
to logg and improve our surface gravity measurement. 
For the handful of WDs in our sample with log g > 7, 
we only fit Balmer lines up to and including HIO since 
the higher order Balmer lines are not observed at higher 
surface gravity. 

The lowest gravity WDs in our sample show Ca and 
Mg lines in their spectra. For reference, the diffusion 
timescale for C a in a T^ff = 10,000 K, 0.2 Mq H-rich 
WD is '^lO^ yr (jPaquette et al.l[l986[ ). Extreme horizon- 
tal branch stars, which can have surface gravities com- 
parable to the lowes t gravity WDs, have longer ^10^ yr 
diffusion timescales (jMichaud et al.l [20081 ). These diffu- 
sion timescales are shorter than the WD evolutionary 
timescale, suggesting there may be on-going accretion 
in these ultra-compact binary systems. Near- and mid- 
infrared observations are needed constrain the possibility 
of accretion. We defer a detailed analysis of the metal 
abundances in ELM WDs to a future paper. For our 
present analysis we exclude the wavelength ranges where 
the metal lines are present in our fits. 

Our error estimates combine the internal error of the 
model fits, obtained from the covariance matrix of the 
fitting algorithm, and the external error, obtained from 
multiple observations of the same object. Uncertainties 
are typically 1.2% in Toff and 0.038 dex in logg (see 
iLiebert et al.ll2005l for details). A measure of the sys- 
tematic uncertainty inherent in the st ellar atmosphere 
models and fitting routines comes from IGianninas et al.l 
(|201lD . who find a systematic uncertainty of « 0.1 dex 
in logg. This is corroborated by the difference we ob- 
serve between the different fitting methods discussed in 
Section 2.1 and here: the mean difference and the dis- 
persion in Toffis 1.1% ± 4.3% and in loggis 0.05 ± 0.12 
dex. Figure [T] plots the properties of the 17 WDs with 
their inter nal error bar s. 

We use iPanei et al.l ()2007[ ) evolutionary tracks to es- 
timate WD mass and luminosity. For purpose of dis- 
cussion, we assume logg < 6 WDs have mass 0.17 Mq 



Brown et al. 



TABLE 1 
White Dwarf Physical Parameters 



Object 



RA 
(h:m:s) 



Dec 
(d:m:s) 



(K) 



logs 

(cm s^^) 



Mass 

(Mrs,) 



90 
(mag) 



Mg 

(mag) 



(kpc) 



J0056-0611 





56:48.232 


-6:11:41.62 


12210 ± 180 


6.167 ±0.044 


0.17 


17.208 ±0.023 


8.0 


0.69 


J0751-0141 


7 


51:41.179 


-1:41:20.90 


15660 ± 240 


5.429 ± 0.046 


0.17 


17.376 ±0.015 


8.0 


0.75 


J0755+4800 


7 


55:19.483 


48:00:34.07 


19890 ± 350 


7.455 ± 0.057 


0.42 


15.878 ±0.019 


9.7 


0.18 


J0802-0955 


8 


02:50.134 


-9:55:49.84 


16910 ± 280 


6.423 ± 0.048 


0.20 


18.604 ±0.012 


8.2 


1.19 


J0811+0225 


8 


11:33.560 


2:25:56.76 


13990 ± 230 


5.794 ± 0.054 


0.17 


18.569 ±0.013 


8.0 


1.30 


J0815+2309 


8 


15:44.242 


23:09:04.92 


21470 ± 340 


5.783 ± 0.046 


0.17 


17.623 ±0.015 


6.7 


1.53 


J0840+1527 


8 


40:37.574 


15:27:04.53 


13810 ± 240 


5.043 ± 0.053 


0.17 


19.141 ±0.018 


8.0 


1.69 


J1046-0153 


10:46:07.875 


-1:53:58.48 


14880 ± 230 


7.370 ± 0.045 


0.37 


17.927 ±0.020 


10.2 


0.36 


J1104+0918 


11:04:36.739 


9:18:22.74 


16710 ± 250 


7.611 ±0.049 


0.46 


16.543 ±0.016 


10.3 


0.18 


J1141+3850 


11:41:55.560 


38:50:03.02 


11620 ±200 


5.307 ±0.054 


0.17 


18.972 ±0.018 


8.0 


1.56 


J1151+5858 


11:51:38.381 


58:58:53.22 


15400 ± 300 


6.092 ±0.057 


0.17 


20.046 ± 0.033 


8.0 


2.57 


J1157+0546 


11:57:34.455 


5:46:45.58 


12100 ±250 


5.054 ±0.071 


0.17 


19.798 ±0.021 


8.0 


2.29 


J1238+1946 


12:38:00.096 


19:46:31.45 


16170 ± 260 


5.275 ±0.051 


0.17 


17.155 ±0.019 


8.0 


0.68 


J1538+0252 


15:38:44.220 


2:52:09.49 


11560 ±220 


5.967 ±0.053 


0.17 


18.528 ±0.015 


8.0 


1.28 


J1557+2823 


15:57:08.483 


28:23:36.02 


12550 ± 200 


7.762 ± 0.046 


0.49 


17.496 ± 0.029 


11.2 


0.18 


J2132+0754 


21:32:28.360 


7:54:28.24 


13700 ± 210 


5.995 ± 0.045 


0.17 


17.904 ±0.019 


8.0 


0.96 


J2338-2052 


23:38:21.505 


-20:52:22.76 


16630 ± 280 


6.869 ±0.050 


0.27 


19.577 ±0.035 


9.0 


1.29 
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Binary Orbital Parameters 











Object 


N„,s 


P 


K 


7 


Spec. Conjunction 


M2 


9 


'''nicrgc 






(days) 


(km s-i) 


(km s-l) 


HJD-2450000 (days) 


(Mq) 




(Gyr) 


J0056-0611 


33 


0.04338 ± 0.00002 


376.9 ±2.4 


4.2 ±1.8 


5864.76305 ± 0.00008 


>0.45 


<0.374 


<0.12 


J0751-0141 


31 


0.08001 ± 0.00279 


432.6 ±2.3 


61.9 ±1.8 


5623.60639 ± 0.00013 


>0.94 


<0.182 


<0.37 


J0755+4800 


26 


0.54627 ±0.00522 


194.5 ±5.5 


42.6 ±3.8 


3730.89702 ± 0.00208 


>0.90 


<0.468 


<28 


J0802-0955 


20 


0.54687 ± 0.00455 


176.5 ±4.5 


27.0 ±3.4 


5623.34011 ± 0.00228 


>0.57 


<0.348 


<79 


J0811+0225 


24 


0.82194 ± 0.00049 


220.7 ±2.5 


77.4 ±1.9 


6329.61559 ± 0.00296 


>1.20 


<0.142 


<160 


J0815+2309 


21 


1.07357 ±0.00018 


131.7 ±2.6 


-37.0 ±2.6 


5623.68831 ± 0.00415 


>0.47 


<0.361 


<620 


J0840+1527 


19 


0.52155 ±0.00474 


84.8 ±3.1 


10.7 ±2.3 


4822.83144 ± 0.00244 


>0.15 


<0.879 


<230 


J1046-0153 


16 


0.39539 ±0.10836 


80.8 ±6.6 


-33.3 ±4.6 


4597.54154 ± 0.00293 


>0.19 


<0.509 


<48 


J1104+0918 


25 


0.55319 ±0.00502 


142.1 ±6.0 


72.5 ±4.0 


3793.84175 ± 0.00269 


>0.55 


<0.831 


<39 


J1141+3850 


17 


0.25958 ± 0.00005 


265.8 ±3.5 


-11.2 ±2.2 


3881.73870 ± 0.00037 


>0.76 


<0.225 


<10 


J1151+5858 


17 


0.66902 ± 0.00070 


175.7 ±5.9 


12.0 ±4.2 


5622.52148 ± 0.00269 


>0.61 


<0.277 


<150 


J1157+0546 


9 


0.56500 ±0.01925 


158.3 ±4.9 


-124.8 ±3.3 


4235.73743 ± 0.00244 


>0.44 


<0.382 


<120 


J1238+1946 


21 


0.22275 ± 0.00009 


258.6 ±2.5 


-6.6 ±1.2 


4236.71749 ± 0.00032 


>0.64 


<0.266 


<7.5 


J1538+0252 


16 


0.41915 ±0.00295 


227.6 ±4.9 


-157.9 ±3.9 


5385.57030 ± 0.00146 


>0.76 


<0.222 


<35 


J1557+2823 


24 


0.40741 ± 0.00294 


131.2 ±4.2 


10.4 ±3.0 


3563.72487 ±0.00127 


>0.43 


<0.886 


<20 


J2132+0754 


35 


0.25056 ± 0.00002 


297.3 ±3.0 


-12.1 ±2.0 


5862.68632 ± 0.00049 


>0.95 


<0.179 


<7.7 


J2338~2052 


25 


0.07644 ± 0.00712 


133.4 ±7.5 


5.2 ±4.8 


5862.75190 ± 0.00049 


>0.15 


<0.554 


<0.95 



Note. — Objects with significant period aliases: J0755-I-4800 (0.349 days), J0840+1527 (0.340 days), J1046-0153 (0.659 days), 
J1104+0918 (0.355 days), J1157+0546 (1.23 days), J1538-I-0252 (0.295 days), and J1557-I-2823 (0.677 and 0.290 days) as seen in 
Figure |3] 



and absolu te magnitude Mg = 8.0 (see lKilic et al.l[2011al : 
iVennes et~al. 2011 ) . We summarize the observed and de- 
rived stellar parameters of the 17 new binaries in Table[TJ 
Positio n and de-reddened g-band magnitude come from 
SDSS (jAihara et al.ll20lTI ) and dheiio is our heliocentric 
distance estimate. 

2.5. Orbital Elements 

We calculate orbital elements and merger times in the 
same way as previous ELM Survey papers, and so we re- 
fer the reader to those papers for the details of our anal- 
ysis. In brief, we measure absolute radi al velocities us- 
the cross-correlation package RVSAO (jKurtz fc MinkI 



[199 



199S) and a high-S/N template. We use the entire spec- 
trum in the cross-correlation. We then use the summed 
spectra (Figure [J) as cross-correlation templates to maxi- 
mize our velocity precision for each individual object. We 
calculate orbital element s by m inimizing x^ for a circular 
orbit using the code of .Kenvon &: Garcia (,1986,1 . Figure 



[3] shows the periodograms for the 17 binaries, and Figure 
0] plots the radial velocities phased to the best-fit orbital 
periods. We use the binary mass function to estimate 
the unseen companion mass; an edge-on orbit with incli- 
nation i = 90° yields the minimum companion mass M2 
and the maximum gravitational wave merger time. 

Table [2] summarizes the binary orbital parameters. 
Columns include orbital period (P) , radial velocity semi- 
amplitude {K), systemic velocity (7), time of spectro- 
scopic conjunction (the time when the object is closest 
to us), minimum secondary mass (M2) assuming i = 90°, 
the maximum mass ratio (g) , and the maximum gravita- 
tional wave merger time Tmcrge- The systemic velocities 
in Table [2] are not corrected for the WDs' gravitational 
redshifts, which should be subtracted from the observed 
velocities to find the true systemic velocities. This cor- 
rection is a few km s~^ for a 0.17 Mq helium WD, com- 
parable to the systemic velocity uncertainty. 
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Fig. 3. — Periodograms for the 17 WD binaries. The best orbital periods have the smallest x values; some binaries are well constrained 
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Fig. 4. — Observed velocities phased to best-fit orbits for the 17 WD binaries (TableO. 



3. RESULTS 

The orbital solutions constrain the nature of the ELM 
WD binaries. Here we discuss the systems with short 
merger times, massive companions, or that may be un- 
der luminous supernovae progenitors. 

3.L 30056-0611 

The ELM WD J0056-0611 has a well-constrained or- 
bital period of 1.0409 ± 0.0005 hr and a semi-amplitude 
of 377 ±2 km s~^. We can calculate its likely companion 
mass if we assume a distribution for the unknown or- 
bital inclination. Although radial velocity detections are 
biased towards edge-on systems (see Section 4), we will 
assume that we are observing a random inclination for 
purpose of discussion. The mean inclination angle for a 
random sample, i = 60°, is then an estimate of the most 
likely companion mass. For J0056— 0611, the most likely 



companion is a 0.61 Mq WD at an orbital separation of 
0.5 i?0. This orbital separation rules out the possibility 
of a main sequence companion. 

There is no evidence for a 0.61 Mq WD in the spec- 
trum of J0056— 0611, but we would not expect there to 
be. If the two WDs in this binary formed at the same 
time, we would expect the 0.61 Mq WD to be 15 - 100 
times less luminous than the 0.17 Mq WD fo r cooling 
ages of 100 Myr - 1 Gyr ([Bergeron et al.lflQQSD . A more 
plausible evolutionary scenario for an ELM WD binary 
like J0056— 0611 is two consecutive phases of common- 
envelope evolution in which the ELM WD is created last, 
giving th e more massive W D a longer time to cool and 
fade fe.g. lKilicet"aI1l2007bD . 

For the most probable companion mass of 0.61 Mq, 
J0056— 0611 will begin mass transfer in 100 Myr. 
IKilic et aLl (|2010bl ) discuss the many possible stellar evo- 
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lution paths for such a system. This system's mass ratio 
o <0.37 suggests t hat mass transfer wih hkely be stable 
(JMarsh et al]|2004D and that J0056-0611 will evolve into 
an AM CVn system. 

3.2. J075 1-0141 

The ELM WD J0751-0141 has a 1.920 ± 0.067 hr or- 
bital period with aliases ranging between 1.85 and 2.05 hr 
(see Figure|3|). The large 433 ±2 km s~^ semi-amplitude 
indicates that the companion is massive, regardless of 
the exact period. The minimum companion mass is 0.94 
Mq. Assuming a random inclination distribution, there 
is a 47% probability that the companion is >1.4 Mq. 

Given that the ELM WD went through a com- 
mon envelope phase of evolution with its companion, 
J0751— 0141 possibly contains a milli-second pulsar bi- 
nary companion. Helium-core WDs are the most com- 
mon t ype of companion i n known milli-second pulsar bi- 
naries (jTauris et al.ll201^ . We have been allocated Cycle 
14 Chandra X-ray Observatory time to search for X-ray 
emission from a possible neutron star. 

For i = 60°, the companion is a 1.32 A/© WD, and 
the system will begin mass transfer in 290 Myr. The ex- 
treme mass ratio q <0.18 means that this system will 
undergo stable mass transfer and evolve into an AM 
CVn system. As the binary orbit widens in the AM 
CVn phase, the mass accretion rate will drop and the 
mass required for the unstable burning of the accreted 
He-layer increases up to several percent of a solar mass. 
The final flash should ignite a thermonucle ar transient 
visible as an underluminou s supernova (Bil dsten et al.l 
l2007HShen fc B ildsten"2009^. It is also possible that the 
helium flash will detonate the massive WD in a double- 
detonation scenario (S im et al.ll2012[ though see Dan et 
al. 2012). If J0751-0141 has a massive WD companion, 
it is a probable supernova progenitor. 

3.3. J0811+0225 

The ELM WD J0811-H0225 has an orbital period of 
19.727 ± 0.012 hr and a semi-amplitude of 221 km s'^ . 
These orbital parameters yield a minimum companion 
mass of 1.20 Mq. That means the companion probably 
exceeds a Chandrasekhar mass. For i — 60°, the most 
likely companion is a 1.70 Mqj neutron star at an orbital 
separation of 4.6 Rq. There is no Fermi gamma-ray de- 
tection at this location, but additional observations are 
needed to determine the nature of this system. 

3.4. J0840+1527 

The ELM WD J0840-H1527 has a best-fit log 5== 
5.043 ± 0.053 near the limit of our model grid. Its best- 
fit orbital period is 12.517 ± 0.114 hr with a significant 
alias at 8.3 hr. Assuming this object is 0.17 Mq, its 
most likely companion is a WD with a comparable mass, 
0.19 Mq, at an orbital separation of 1.9 Rq. If, on the 
other hand, J0840-I-1527 were a 3 Mq main sequence 
star, its companion would have an orbital separation of 
4.3 Rq - a separation comparable to the radius of a 3 
Mq star, and thus physically implausible. There is no 
evidence for mass transfer in this system. We conclude 
that J0840+1527 is a pair of ELM WDs. 

3.5. J1141+3850, J1157+0546, and J1238+1946 



J1141-)-3850, J1157-I-0546, and J1238-H1946 are the 
other systems containing WDs near the low-gravity limit 
of our stellar atmosphere model grid, but their k = 158 - 
266 km s~^ semi-amplitudes are significantly larger than 
that of J0840-I-1527. If we assume that the objects are 
0.17 Mq ELM WDs, then the binary companions have 
minimum masses of 0.45 - 0.75 Mq. If we instead as- 
sume that the objects are main sequence stars, then the 
required orbital separations are comparable to the ra- 
dius of the main sequence star and physically impossible. 
There is no evidence for mass transfer in these systems. 

We conclude that J1141+3850, J1157-H0546, and 
J1238-hl946 are ELM WDs with likely WD companions. 
For J1141-I-3850, there is a 36% probability that the com- 
panion is >1.4 Mq, possibly a milli-second pulsar. For 
J1141-I-3850 and J1238-I-1946, mass transfer wiU begin in 
7-10 Gyr, making them AM CVn progenitors and possi- 
ble underluminous supernovae progenitors. 

3.6. J2132+0754 

The ELM WD J2132-H0754 has a well-constrained or- 
bital period of 6.0134 ± 0.0004 hr and semi-amplitude of 
297 ± 3 km s^^. The minimum companion mass is 0.95 
Mq, and there is a 48% probability that the companion 
is >1.4 Mq, possibly a milli-second pulsar. For i = 60°, 
the most likely companion is a 1.33 Mq WD that will 
begin mass transfer in 6 Gyr. That makes J2132-I-0754 
a likely AM CVn progenitor and another possible under- 
luminous supernovae progenitor. 

3.7. J 2338- 2052 

The ELM WD J2338-2052 has an orbital period of 
1.834 ± 0.170 hr and a semi- amplitude of 133 ± 7 km 
s^^. In this case the companion is another ELM WD; 
for i — 60°, the most likely companion is a 0.17 Mq WD 
at an orbital separation of 0.56 Rq. Given the unity 
mass ratio, this system will undergo unstable mass trans- 
fer and will merge to form a single ~0.4 Mq WD. This 
system will merge in less than 1 Gyr. 

4. DISCUSSION 

With these 17 new discoveri es, plus the pu l sating ELM 
WD discovery published by iHermes et al.l (|2013[ ). the 
ELM Survey has found 54 detached, double degenerate 
binaries; 31 of the binaries will merge within a Hubble 
time. Table [3] summarizes the properties of the systems. 
Eighty percent of the ELM Survey binaries are formally 
ELM WD systems. 

4.1. Significance of Binary Detections 

We find low mass WDs in compact binaries, binaries 
that must have gone through common envelope evolu- 
tion. This makes sense because extremely low mass WDs 
require significant mass loss to form; the Universe is 
not old enough to produce an extremely low mass WD 
through single star evolution. Yet four objects published 
in the ELM Survey have no significant velocity variabil- 
ity, two of which are ELM WDs (see Table [3]). To under- 
stand whether or not these stars are single requires that 
we understand the significance of our binary detections. 

Each ELM Survey binary is typically constrained by 
10-30 irregularly spaced velocities with modest errors. 
Given that we determine orbital parameters by minimiz- 
ing x^, the i^-test is a natural choice. We use the F-test 
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Fig. 5. — Orbital period vs. semi-amplitude for the ELM Survey 
binaries. Previously published binaries are drawn with solid red 
circles; the 17 binaries from this paper are drawn with solid blue 
diamonds. Dashed lines indicate approximate companion mass un- 
der the assumption that A/i = 0.2 Mq and i = 60°. 

to check whether the variance of the data around the or- 
bital fit is consistent with the variance of the data around 
a constant velocity (we use the weighted mean of the ob- 
servations), i^-test probabilities for the published ELM 
Survey binaries are < 0.01. In other words, our binaries 
have significant velocity variability at the > 99% confi- 
dence level. 

In the null cases we need to calculate the likelihood of 
not detecting a binary. This is a trickier problem, and 
one that we approach with a Monte Carlo calculation. 
We start by selecting a set of observations (times, ve- 
locity errors) and an orbital period and semi- amplitude. 
We convert observation times to orbital phases using a 
randomly drawn zero time, and calculate velocities at 
those phases summed with a randomly drawn velocity 
error. We perform the F-test, using 0.01 as a detection 
threshold. We repeat this calculation 10,000 times for a 
given orbital period and semi- amplitude, and then select 
a new orbital period and semi-amplitude to iterate on. 
This analysis is done for each object. 

We find that the datasets for our 17 new binaries have 
a median 99.9% likelihood of detecting K — 200 km s^^ 
binaries, a 97% likelihood of detecting K — 100 km s~^ 
binaries, and a 44% likelihood of detecting K = 50 km 
s^^ binaries. It is no surprise that we are less likely 
to detect a low semi-amplitude binary, but this analysis 
suggests that we can be quite confident of detecting K > 
100 km s~^ systems. We find very similar likelihoods for 
detecting binaries containing ELM WDs in the full ELM 
Survey sample. 

The datasets for the null cases typically contain fewer 
observations and so are not as well-constrained. For 
J0900+0234, a .16 M^t, ELM WD with no observed ve- 
locity variation (| Brown et a l. 2012b), the likelihoods of 
detecting a. K ^ 200, 100, and 50 km s^^ binary are 87%, 
57%, and 10%, respectively. Additional observations are 
required to claim this ELM WD as non- variable. 

There is, of course, an orbital period dependence to the 
detections, and periods near 24 hr are the most prob- 













Ml| 








2 


• ELM 


Su 


rvey 












"♦ This 


paper 




























♦ " 










J1741 








J0611 _ 










• 








J0755 


















Jai32 * 


h 


~ Mchan 














♦ - 




ra 






















• 




• 


• 


_ 










• 








%♦ ♦ 
















* * 


1 


• 






• 

• • 

• ♦ 




• 


• 


• 

• * 


5 


- J065i 




• 


J0056 






• 










♦ 
J2338 


• 


♦ " 


















• 




" 




Mill 


1 1 1 1 1 


III 






1 1 1 



10'* 



IC 



10" 



10» 

[yr] 



10'° 10" 



Fig. 6. — Gravitational wave merger time versus total system 
mass for the ELM Survey. We calculate system mass assuming 

o 

j = 60 when orbital inclination is unknown. Previously published 
ELM Survey binaries are drawn with solid red circles, and the 17 
new binaries from this paper are drawn with solid blue diamonds. 
Six binaries have probable masses exceeding the Chandrasekhar 
mass; three have merger times less than 10 Gyr, as indicated by 
the dotted line. The minimum companion mass for J0811+0225 is 
1.2 Mq. 

lematic. Taken together, our datasets have a median 
39% likelihood of detecting a K — 100 km s^^ binary at 
P = 24 hr. Yet we remain sensitive to longer periods: 
our datasets have median 99% and 98% likelihoods of 
detecting a K = 100 km s~^ binary at P = 18 hr and 
P = 36 hr, respectively. 

Figure [5] plots the observed distribution of P and K 
for the ELM Survey binaries. The dashed lines indicate 
the approximate companion mass assuming Mi is 0.2 Mq 
and i is 60°. At K = 100 km s^^, we can detect compan- 
ion masses down to 0.1 Mq at 2 hr orbital periods and 
0.55 Mq companions at 2 day orbital periods. Systems 
with K < 100 km s~^ are the realm of ^0.2 Mq com- 
panions, and we observe a half-dozen systems with these 
parameters. Our incompleteness at K < 100 km s~^ im- 
plies there are quite likely more ELM WDs with ~0.2 Mq 
companions; the remaining ELM WD candidates that do 
not show obvious velocity change in a couple observations 
are possible examples of such low amplitude systems. 

Finally, orbital inclination acts to increase the diffi- 
culty of identifying lone ELM WDs. Consider the set of 
46 ELM WDs in Table [3] with logg < 7, two of which 
are non- variable. Their median K is 240 km s~^, a semi- 
amplitude that would appear < 50 km s^^ at i < 12°. 
If the 46 objects have randomly distributed orbital in- 
clinations, one object should have i < 12° and a second 
i < 17°. We conclude there is no good evidence for a 
lone ELM WD in our present sample. This is in stark 
contrast to the population of 0.4 Mq WD s in the solar 
neighb orhood, of which 20%-30% are single (jBrown et al.l 
l2011aD . 

4.2. The Future of ELM WDs 

One of the most exciting aspects of our ELM WD bi- 
naries is that many have gravitational wave merger times 
less than a Hubble time, and in one case as short as 1 
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Fig. 7. — Gravitational wave merger time versus mass ratio 

q for the ELM Survey. We calculate mass ratio assuming i = 

o 
60 when orbital inclination is unknown. Previously published 

ELM Survey binaries are drawn with solid red circles, and the 17 

new binaries from this paper are drawn with solid blue diamonds. 

Dotted lines mark the approximate thresholds for stability of mass 

transfer during the Roche lobe overflow phase (Marsh ot al. 2004), 

drawn for systems that will merge in <10 Gyr. Stable mass transfer 

systems will evolve into AM CVn systems; unstable mass transfer 

systems will possibly merge. 



Myr (jBrown et al.ll2011cl: iHermes et al.ll2012bf ). A nat- 
ural question, then, is what will happen when the WDs 
merge. Notably, 5 of our 17 new binaries have probable 
system masses (for i = 60°) in excess of a Chandrasekhar 
mass. 

Figure [6] plots the distribution of merger time and sys- 
tem mass for the ELM Survey binaries assuming i = 60° , 
except when inclination is known. The majority of our 
binaries have probable masses below the Chandrasekhar 
mass and are not potential supernova Type la progeni- 
tors, but six systems appear to have either neutron star 
or massive WD companions. The evolution of these six 
systems depends on the stability of mass transfer dur- 
ing the Roche lobe overflow phase, and thus on the mass 
ratio of the binary. 

Figure[7]plots the distribution of merger time and mass 
ratio q for the ELM Survey binaries, again assuming 
i = 60° except when inclination is known. The Chan- 
drasekhar mass systems generally have extreme q < 0.15 
and so will evolve into stable mass transfer AM CVn sys- 
tems. If the accretors are massive WDs, these systems 
should experience a large helium flash that may appear 
as an imderluminous .la sup ernova (e.g. Bildstcn et al. 
I2007t iShen fc Bildstenl |2009[ ). It is also possible that 
the helium flash will deto nate the massiv e WD in a 
double-detonation scenario ()Sim et al.ll2012f) . The Chan- 
drasekhar mass systems are thus potential supernovae 
progenitors. However, the outcome of a me rger of an 
ELM WD with a massive WD is uncertain (see lDan et alJ 
[20T1 . 

Systems with q > 2/3, such as J2338— 2052, will expe- 
rience unstable m ass transfer. These systems may merge 
(|Dan et al.ll2011[ ) to form a single low-mass WD, R CrB 
star, or helium-burning sdB star. Systems with interme- 
diate q ~ 1/2, such as J0651+2844, will experience either 



stable or unstable mass transfer depending on their spin- 
orbit coupling. After we constrain the remaining ELM 
WD candidates plotted in Fig. [1] we look forward to cal- 
culating their space density and merger rate. ELM WDs 
binaries are clearly an important channel for AM CVn 
star formation, and thus an important source of strong 
gravitational waves in the mHz regime. 

5. CONCLUSION 

We perform stellar atmosphere fits to the entire col- 
lection of single-epoch spectra in the HVS Survey and 
identify 57 new low mass WD candidates. Follow-up 
spectroscopy reveals 17 WDs with significant velocity 
variability, 12 of which are ELM WDs. Presently, the 
ELM Survey sample is consistent with all of the ELM 
WDs being part of close binaries. ELM WDs are thus 
signposts for binaries that are strong gravitational wave 
sources and possible supernovae progenitors. 

Four binaries in our sample contain logg ~ 5 objects, 
objects that might not be considered WDs or proto-WDs 
if not for their observed orbital motion. The existence 
of logg ~ 5 WDs motivates us to de velop a new set 
of ste llar atmosphere models following iGianninas et alJ 
(j2011h . Interestingly, all of the lowest surface gravity 
(logg < 6) WDs in our sample have metal lines in their 
spectra. This issue will be studied in detail in a future 
paper; infrared observations are needed to constrain the 
possibility of debris disks as the source of Ca and Mg 
accretion in these compact binary systems. 

Recent disc overies of pulsations in sever al of our ELM 
WD targets (jHermes et all l2012d I2013D enable us to 
probe the interiors of low-mass, presumably He-core 
WDs using the tools of asteroseismology. Due to the 
HVS color selection, all of the newly identified systems 
are hotter than 10,000 K and therefore too hot to show 
the g-mode pulsati ons detected in cooler ELM WDs 
(|C6rsico et al.llMl) . 

Four binaries in our sample have massive ^0.9 Mq 
companions and <10 Gyr merger times. If the unseen 
companions are massive WDs, these extreme mass ratio 
binaries will undergo stable mass transfer and will evolve 
into AM CVn systems and potentially future .la or un- 
derluminous supernovae. Thus the systems J0751— 0141, 
J114H-3850, J1238+1946, and J2132+0754 are possible 
supernova progenitors. 

The other possibility is that the massive compan- 
ions in these four binaries are neutron stars. Given 
their past common envelope evolution, the systems quite 
possibly contain milli-second pulsars. A fifth system, 
J081 1-1-0225, appears to have a minimum companion 
mass close to a Chandrasekhar mass. We are conducting 
follow-up optical, radio, and X-ray observations to es- 
tablish the nature of these extreme mass ratio ELM WD 
binaries. Optically visible WD companions of neutron 
stars are useful for constraining the binary mass ratio 
and also to calibrate the spin-down ages of milli-second 
pulsars. 

Given our present observations, we expect that the 
ELM Survey will grow to 100 binaries in a few years. A 
large sample is important for finding more systems like 
J0651-I-2844, a 12.75 min period system that provides us 
with a laboratory for measuring the spin-orbit coupling 
and tidal heating of a rapidly merging pair of WDs. We 
expect that photometric observations will discover more 
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pulsating ELM WDs. Finally, a larger sample of ELM 
WDs will be an important source of gravitational wave 
verification sources for eLISA and other future gravita- 
tional wave detection experiments. 

We thank M. Alegria, P. Canton, S. GotiUa, J. McAfee, 
E. Martin, A. Milone, and R. Ortiz for their assistance 
with observations obtained at the MMT Observatory, 
and P. Berlind and M. Calkins for their assistance with 
observations obtained at the Fred Lawrence Whipple Ob- 
servatory. We also thank P. Bergeron and D. Koester 
for their invaluable assistance in the computation of our 
stellar atmosphere model grids. This research makes use 
of the SAO/NASA Astrophysics Data System Biblio- 
graphic Service. This project makes use of data prod- 
ucts from the Sloan Digital Sky Survey, which has been 
funded by the Alfred P. Sloan Foundation, the Partic- 
ipating Institutions, the National Science Foundation, 
and the U.S. Department of Energy Office of Science. 
This work was supported in part by the Smithsonian In- 
stitution. 

Facilities: MMT (Blue Channel Spectrograph), 
FLWO :1.5m (FAST) 



10 



Brown et al. 



TABLE 3 

Merger and Non-Merger Systems in the ELM Survey 



Object 


TcB 


logg 


P 


K 


Mass 


Ma 


A/a (60°) 


■^mcrgc 


Ref 




(K) 


(cm s ^) 


(days) 


km s^^ 


Mq 


Mq 


Mq 


Gyr 




J0022-1014 


18980 


7.15 


0.07989 


145.6 


0.33 


>0.19 


0.23 


<0.73 


6 


J0056-0611 


12210 


6.17 


0.04338 


376.9 


0.17 


>0.46 


0.61 


<0.12 




J0106-1000 


16490 


6.01 


0.02715 


395.2 


0.17 


0.43 




0.037 


7 


J0112+1835 


9690 


5.63 


0.14698 


295.3 


0.16 


>0.62 


0.85 


<2.7 


1 


J0651+2844 


16530 


6.76 


0.00886 


616.9 


0.26 


0.50 




0.0011 


3,15 


J0751-0141 


15660 


5.43 


0.08001 


432.6 


0.17 


>0.94 


1.32 


<0.37 




J0755+4906 


13160 


5.84 


0.06302 


438.0 


0.17 


>0.81 


1.12 


< 0.22 


2 


J0818+3536 


10620 


5.69 


0.18315 


170.0 


0.17 


>0.26 


0.33 


<8.9 


2 


J0822+2753 


8880 


6.44 


0.24400 


271.1 


0.17 


>0.76 


1.05 


<8.4 


4 


J0825+1152 


24830 


6.61 


0.05819 


319.4 


0.26 


>0.47 


0.61 


<0.18 





J0849+0445 


10290 


6.23 


0.07870 


366.9 


0.17 


>0.64 


0.88 


<0.47 


4 


J0923+3028 


18350 


6.63 


0.04495 


296.0 


0.23 


>0.34 


0.44 


<0.13 


2 


J1005+0542 


15740 


7.25 


0.30560 


208.9 


0.34 


>0.66 


0.86 


<9.0 





J1005+3550 


10010 


5.82 


0.17652 


143.0 


0.17 


>0.19 


0.24 


< 10.3 





J1053+5200 


15180 


6.55 


0.04256 


264.0 


0.20 


>0.26 


0.33 


<0.16 


4,9 


J1056-I-6536 


20470 


7.13 


0.04351 


267.5 


0.34 


>0.34 


0.43 


< 0.085 





J1112+1117 


9590 


6.36 


0.17248 


116.2 


0.17 


>0.14 


0.17 


< 12.7 


16 


J1141+3850 


11620 


5.31 


0.25958 


265.8 


0.17 


>0.76 


1.05 


<9.96 




J1233+1602 


10920 


5.12 


0.15090 


336.0 


0.17 


>0.86 


1.20 


<2.1 


2 


J1234-0228 


18000 


6.64 


0.09143 


94.0 


0.23 


>0.09 


0.11 


< 2.7 


6 


J1238+1946 


16170 


5.28 


0.22275 


258.6 


0.17 


>0.64 


0.88 


< 7.49 




J1436+5010 


16550 


6.69 


0.04580 


347.4 


0.24 


>0.46 


0.60 


<0.10 


4,9 


J 1443+ 1509 


8810 


6.32 


0.19053 


306.7 


0.17 


>0.83 


1.15 


<4.1 


1 


J1630+4233 


14670 


7.05 


0.02766 


295.9 


0.30 


>0.30 


0.37 


< 0.031 


8 


J1741+6526 


9790 


5.19 


0.06111 


508.0 


0.16 


> 1.10 


1.55 


<0.17 


1 


J1840+6423 


9140 


6.16 


0.19130 


272.0 


0.17 


>0.64 


0.88 


<5.0 


1 


J2103-0027 


10000 


5.49 


0.20308 


281.0 


0.17 


> 0.71 


0.99 


<5.4 





J2119-0018 


10360 


5.36 


0.08677 


383.0 


0.17 


>0.75 


1.04 


<0.54 


2 


J2132+0754 


13700 


6.00 


0.25056 


297.3 


0.17 


>0.95 


1.33 


< " 


'.70 




J2338-2052 


16630 


6.87 


0.07644 


133.4 


0.27 


>0.15 


0.18 


<C 


.95 




NLTT11748 


8690 


6.54 


0.23503 


273.4 


0.18 


0.76 




7 


2 


5,10,11 


J0022+0031 


17890 


7.38 


0.49135 


80.8 


0.38 


> 0.21 


0.26 






6 


J0152+0749 


10840 


5.80 


0.32288 


217.0 


0.17 


>0.57 


0.78 






1 


J0730+1703 


11080 


6.36 


0.69770 


122.8 


0.17 


>0.32 


0.41 









J0755+4800 


19890 


7.46 


0.54627 


194.5 


0.42 


>0.90 


1.18 








J0802-0955 


16910 


6.42 


0.54687 


176.5 


0.20 


>0.57 


0.76 








J0811+0225 


13990 


5.79 


0.82194 


220.7 


0.17 


> 1.20 


1.71 








J0815+2309 


21470 


5.78 


1.07357 


131.7 


0.17 


>0.47 


0.63 








J0840+1527 


13810 


5.04 


0.52155 


84.8 


0.17 


>0.15 


0.19 








J0845+1624 


17750 


7.42 


0.75599 


62.2 


0.40 


>0.19 


0.22 









J0900+0234 


8220 


5.78 




< 24 


0.16 










1 


J0917+4638 


11850 


5.55 


0.31642 


148.8 


0.17 


>0.28 


0.36 






12 


J1046-0153 


14880 


7.37 


0.39539 


80.8 


0.37 


>0.19 


0.23 








J1104+0918 


16710 


7.61 


0.55319 


142.1 


0.46 


>0.55 


0.70 








J1151+5858 


15400 


6.09 


0.66902 


175.7 


0.17 


>0.61 


0.84 








J1157+0546 


12100 


5.05 


0.56500 


158.3 


0.17 


>0.45 


0.59 








J1422+4352 


12690 


5.91 


0.37930 


176.0 


0.17 


> 0.41 


0.55 






2 


J1439+1002 


14340 


6.20 


0.43741 


174.0 


0.18 


>0.46 


0.62 






2 


J1448+1342 


12580 


6.91 




<35 


0.25 










2 


J1512+2615 


12130 


6.62 


0.59999 


115.0 


0.20 


>0.28 


0.36 






2 


J1518+0658 


9810 


6.66 


0.60935 


172.0 


0.20 


>0.58 


0.78 






1 


J1538+0252 


11560 


5.97 


0.41915 


227.6 


0.17 


>0.77 


1.06 








J1557+2823 


12550 


7.76 


0.40741 


131.2 


0.49 


>0.43 


0.54 








J1625+3632 


23570 


6.12 


0.23238 


58.4 


0.20 


>0.07 


0.08 






6 


J1630+2712 


11200 


5.95 


0.27646 


218.0 


0.17 


>0.52 


0.70 






2 


J2252-0056 


19450 


7.00 




<25 


0.31 










2 


J2345-0102 


33130 


7.20 




<43 


0.42 










2 


LP400-22 


11170 


6.35 


1.01016 


119.9 


0.19 


> 0.41 


0.52 
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APPENDIX 
DATA TABLE 



Table [4] presents our radial velocity measurements. The Table columns include object name, heliocentric Julian date 
(based on UTC), heliocentric radial velocity (uncorrected for the WD gravitational redshift), and velocity error. 
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TABLE 4 
Radial Velocity Measurements 



Object 


HJD 


^helio 




(days-2450000) 


(km s-1) 


J0056-0611 


5864.778351 


308.7 ±9.7 




5864.786174 


-64.5 ±9.5 




5864.788281 


-162.4 ±10.7 




5864.789600 


-249.1 ± 8.8 




5864.790804 


-258.9 ±10.8 




5864.792007 


-304.9 ± 10.3 



Note. — (This table is available in its en- 
tirety in machine-readable and Virtual Observa- 
tory forms in the online journal. A portion is 
shown here for guidance regarding its form and 
content.) 



